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Abstract. We present experimental observations of diffusion resonances for the 
quantum kicked rotor with weak decoherence. Cold caesium atoms are subject to 
a pulsed standing wave of near-resonant light, with spontaneous emission providing 
environmental coupling. The mean energy as a function of the pulse period 
is determined during the late-time diffusion regime for a constant probability of 
spontaneous emission per unit time. Structure in the late-time energy is seen 
to increase with physical kicking strength. The observed structure is related to 
Shepelyansky's predictions of the initial quantum diffusion rates. 
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1. Introduction 

The atom-optics realization of the kicked rotor has enabled the experimental study of 
the transition between quantum and classical behaviour for this fundamental non-linear 
system. For example, the effects of decoherence, the mechanism whereby quantum 
interference effects are destroyed via environmental coupling PJ, have been studied in 
the quantum system. More classical-like behaviour is observed when decoherence is 
added, either by spontaneous emission events [21 HUE!, amplitude noise IE] or timing 
noise [7J [H] . However, recent theoretical studies have concentrated on what is perhaps 
a more direct approach to studying the quantum-to-classical transition - varying the 
action of the system and thereby the effective Planck's constant, i.e., increasing the 
action to investigate the limit in which 'h'— > [3 EH|. Of particular interest is the 
behaviour found for intermediate values of 'h', for which regions of enhanced diffusion 
(known as diffusion "resonances") are predicted to exist by simulations and analytical 
results. 

These theoretical studies have concentrated on the case where the classical 
stochasticity parameter of the system stays the same and only the effective Planck's 
constant is varied. However, in cold atom kicked rotor experiments, the most accessible 
parameter regime exists for the situation where the power in the kicking laser is held 
constant instead, while the pulsing period (which is proportional to the effective Planck's 
constant) is changed. Numerically, it is found that this gives rise to different but 
analagous diffusion resonances to those found in the aforementioned simulations. Here, 
we present observations of diffusion resonances in this experimentally accessible regime, 
the structure of which can be traced to a scaling formula for the initial quantum 
diffusion rate derived by Shepelyansky jTT]. The experimental results exhibit markedly 
different behaviour to that predicted by classical calculations. An increase in the 
complexity of the resonance structure is seen as the physical kicking strength is increased. 
Experimentally, the resonances are slightly diminished due to the non-uniform instensity 
profile of the kicking laser but the expected structure is still clearly visible. 

The resonances we observe are of a different nature to the quantum resonances 
previously studied in the atom optics kicked rotor f, showing non-trivial dependence on 
the kicking strength as well as the scaled Planck's constant. In particular, we note that 
the emphasis of our study is different from that of the related work by d'Arcy et al. 
which focused on such quantum resonance behaviour and not on the diffusion resonance 
behaviour at intermediate values of h. Additionally, we note the difference between our 
experiments and those of Klappauf et al. p2] in which anomalous diffusion behaviour 
was studied as a function of kicking strength for just a few values of the pulsing period. 
Our investigation is essentially the converse of this experiment: Energies are measured 
as a function of the pulsing period for just a few values of a kicking strength parameter, 
and display a resonance structure which cannot be inferred from previous results. 

| The simplest quantum resonances (which are also those that have been previously studied 
experimentally) occur when the scaled Planck's constant is an integer multiple of 2ir. 
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The system and model that we study is presented in Sec. II, while the analysis of 
classical and quantum momentum diffusion in this system is discussed in Sec. III. 
Decoherence effects due to the small rates of spontaneous emission present in our 
experiments are also considered here. Our experimental set-up is described in Sec. 
IV, and the results (experimental, theoretical quantum and classical calculations) are 
presented in Sec. V. The conclusion is in Sec. VI. 

2. System and Model 

For our system we use a laser-cooled cloud of caesium atoms of initial temperature ~ 20 

/j,K interacting with a standing wave of off-resonant laser light. The laser is pulsed with 

period T and pulse profile f(t). If the detuning of the laser from the atomic transition 

is large, the internal atomic dynamics can be eliminated and the motion of the caesium 

atoms is described by the single particle Hamiltonian jHj 

•f? hO N 
H= ^-^cos(2k L x) Y,f(t~nT), (I) 

2m 8 n=0 

where x and p are operators representing the atomic position and momentum, 
respectively, and is the wave number of the laser light. The effective potential 
strength, Q c s = Vt 2 {s^/8^ + s^/5^ + 543/^43), accounts for dipole transitions between 
different combinations of hyperfine levels in the caesium atoms (6Si/2(-F = 4) — > 
6P3/2(i r " = 5)), where 5ij are the corresponding detunings between laser and atomic 
transiton frequencies, and Q/2 is the resonant single-beam Rabi frequency. If we assume 
equal populations of atoms in all ground state Zeeman sublevels, then s 45 = s 44 = ^, 
and S43 = ^g. It is useful to rewrite this Hamiltonian in appropriate dimensionless units 
as 

-2 N 

W = £--*cos$X;/(T-n) l (2) 

1 n=0 

which is the usual expression for the Hamiltonian of the standard kicked rotor system. 
In these units - which will be referred to as "scaled units" - the position operator is 
defined by <p = 2k^x, the momentum operator is p = 2k^Tp/m, time is rescaled as 
t = t/T, and our new Hamiltonian is related to Eq. (0} by H = (AklT 2 /m)H. The 
classical stochasticity parameter is given by k = fl e su r TT p , where r p is the pulse length 
in unsealed time and u r = hk1/2m. In our experiments, /(r) is a good approximation 
to a square pulse, i.e., /(r) = 1 for < r < a, where a = t p /T. Note that k = tz/a. 

In scaled units we have [4>, p] = ik, with k = 8u r T, so that the quantum behaviour 
of our system is reflected by an effective Planck's constant, k , which increases as we 
increase the pulse period T. This reflects our ability to change the total action in the 
system, and hence how classically our system behaves (for larger k values the quantum 
nature of the system should be more apparent). Note that the effective Planck's constant 
is proportional to the ratio of the total classical action of the system to h. 

The natural experimental unit for momentum is that of two photon recoils, 2hk^, 
and pj (2hkj J ) will henceforth be referred to as the momentum in experimental units. We 
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note the relationship p/k = p/(2hk L ) and also define the quantity 0d = K /k = Q eS r p /S 
as a dimensionless measure of the physical kicking strength. Experimentally, it is easier 
to hold this quantity constant, rather than k, as T is varied, because a constant value of 
0d corresponds to constant pulse duration, standing wave detuning and power (whereas 
k is proportional to T). 

Our system is coupled to its environment via atomic spontaneous emission events, 
which occur when the caesium atoms absorb photons from the standing wave [2] and 
then spontaneously re-emit the photons in random directions. We characterise the level 
of this decoherence by the probability of spontaneous emission per atom per kick, 77. 
Given the large detuning used in our experiments, i.e., Cleg/ 5 <C 1, this process may be 
modeled by the following master equation for the density operator w of the system [3], 

i - n N 

w = - -[H, w]--£/(r- n)[cos 2 (0/2), w] + 

* « n=0 

n N fl 
+ 2^£/(r-7i) / duN{u)e iu(t>/2 

a n=0 

x cos(0/2)u;cos(0/2)e- i ^ /2 , (3) 

where N(u) is the distribution of recoil momenta projected onto the axis of the standing 
wave, and [.,.]+ denotes an anti-commutator. Simulations of Eq. Q are used for 
comparisons with the experiment. 

3. Momentum Diffusion 

We measure the total kinetic energy of the cloud after N kicks, which depends on the 
initial energy of the cloud plus the increase in the kinetic energy resulting from the 
kicks. The amount of increase for kick number n is the momentum diffusion rate, given 
by 2D{n) = (Pn +1 ) — (p 2 ), where we denote p = p(t' = 0), p\ = pit' = 1), etc.. For a 
kicked rotor system with a sufficiently broad initial momentum distribution, we expect 
D(0) = D(l) = k 2 /4. The system then passes through an initial quantum diffusion 
period lasting typically for around 5 kicks jTUj, with a diffusion rate approximated by 
the result of Shepelyansky (under the conditions k > 1 and k ^> k) [TTj . 

D, = ^--J 2 (K q )-j' 1 (K^ + Jl(K q ) + J!(K cl )j , (4) 

where K q = 2Ksm{k/2)/k. Note that the long time classical diffusion rate is also given 
by Eq. (pjj, but with K q -> k (i.e., k ->0). 

Without decoherence, the system generally settles into a localised state [TH], but 
the loss of phase coherence produced by the addition of spontaneous emission causes 
the system to settle instead into a final steady state diffusion regime, with a late time 
diffusion rate which may be approximated by the formula [3 EH Ho] 

00 

A* = I>(l-*7) n A>(n), (5) 

n=0 
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where D (n) is the diffusion rate at the nth kick for a kicked rotor without decoherence. 
Essentially, this formula assumes that dynamical correlations over particular time 
intervals which give rise to the late time diffusion rates are suppressed by a factor 
equal to the probability that a spontaneous emission occurs within that time interval. 
The correlations taken over a set number of kicks give rise to the diffusion rates seen 
in the kicked rotor without decoherence after that number of kicks, which leads to the 
late time diffusion rate being an appropriate weighted average over the diffusion rates 
as the kicked rotor "settles down" Thus, the diffusion rates in the first few kicks 
are essentially "locked in" by the spontaneous emission events, and we observe similar 
structure in the late time diffusion rates as we vary T to that observed in the initial 
quantum diffusion rates. 

The structure we observe in the initial quantum diffusion rates as we vary T 
for constant <pd with diffusion rates measured in experimental units is particularly 
interesting. We can express Shepelyansky's formula in this regime as 

with K' q = 2<pd sin(4o; r T). We then see that any structure in the diffusion rates is 
periodic in T with period 27r/8o; r (In fact, from Eq. (j5J) this is also true for the late time 
diffusion rates). We also see that the form of the structure depends solely on the value 
of </><j. Fig. n shows the initial quantum diffusion rate as a function of pulse period for 
varying values of 0d- We see the regular feature of a peak at the quantum resonance 
when k = 2tt (T = 60.4 us), and we see increasing numbers of pronounced peaks or 
diffusion resonances as we increase the value of 0d- 

The classical diffusion rate can be similarly found in this regime to be 

D'^^i—H^-Jl^ + Jj^ + Jl^ , (7) 



2 V2 

with k = 8u) T T<j)d- The dotted curves in Fig. Q show the classical rates for various 
values of 0d- These rates oscillate around the quasilinear value which in these units is 
(0d) 2 /4, with the oscillations increasing in frequency with increasing 0d- For any given 
4>d however, the structure in the classical diffusion rate is markedly different to that in 
the initial quantum diffusion rate. 

In the experiment, late-time energies were measured and not initial diffusion rates. 
However, the diffusion resonances discussed above are still observable. The existence 
of the same structures in the late time diffusion rates (in the presence of decoherence) 
as those in the initial quantum rates has been verified using the simulations described 
later in this paper. An approximate formula for the energy after N kicks can be found 
by summing the diffusion rate at each kick, that is 

((p^I) 2 ) N-l 

E'(N) = VmiJ = J2 D '( n , T )- (8) 

2 n =0 

Therefore, with the initial and final diffusion rates both displaying these structures, it 
follows that the energy at the iVth kick should also display them. 
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Figure 1. Theoretical initial momentum diffusion rates in the quantum case (solid 
line), D' q as given by Eq. JSJl, and the classical case (dotted line), D' cl in Eq. 10. The 
rates are given in experimental units as a function of T for (a) <pd — 3, (b) <pd = 4.5, 
(c) <fid = 6 and (d) 4>d — 7.5. The quantum resonance is seen in the quantum diffusion 
rate at T = 60.4 us. 



4. Experimental Setup 

The experimental setup used was much the same as that used previously in our 
quantum chaos investigations ^21 E] with a few modifications. A standard six-beam 
magneto-optical trap (MOT) was used to trap and cool approximately 10 5 caesium 
atoms. The trapping laser frequency was set about 10 MHz to the red of the 
6Si/2(-F = 4) — > 6P 3 / 2 (F' = 5) transition. A second (repump) laser was locked to 
the 6Si/2(-F = 3) — ► GP^^F' = 4) transition to return those atoms lost to the F = 3 
ground state to the trapping cycle. After a 20 ms cooling phase prior to kicking, the 
cloud had a temperature of approximately 20 fiK and a width of o" c i ou d ~ 270 //m in 
its position distribution. Kicking of the cloud occurred for up to 2 ms during the 10 
ms free expansion phase, at the completion of which the cloud was "frozen" in space by 
the molasses beams and imaged. The repumping beam was left on during the kicking 
to prevent loss of atoms to the F = 3 ground state. The resultant heating effect was 
negligible for our experiments. 

A third laser was used to create a pulsed optical standing wave across the cloud. 
A relatively high power laser diode was injection locked with a frequency stabilised 
external cavity laser, giving a beam of up to 22 mW CW power. For fast switching the 
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beam passed through a 80 MHz Acousto-Optic Modulator (AOM) in front of a single 
mode polarisation preserving optical fibre. Temporal modulation was provided via the 
RF supply to the AOM, generating pulse shapes very close to rectangular. The linearly 
polarised beam was then collimated resulting in a beam radius at the cloud of 2o"b e am = 1 
mm. Finally, to create a standing wave the beam was retroreflected by a mirror outside 
the vacuum cell. The atoms experienced a range of optical potential depths as the cloud's 
width was comparable in size to that of the laser beam. If 0d,max is the kicking strength 
along the beam axis then the mean value was found to be 0d,mean ~ O.770d,max with a 
standard deviation of 18%. This effect and additional factors leading to a slightly greater 
spread in 4> Tm d values is discussed below. In the following, 0d will always refer to 0d,mean- 
The detuning of the kicking beam to the blue of the F = 4 — > F' = 5 transition was 
monitored by overlapping the trapping and kicking beams and monitoring the resultant 
beat signal. Both the beam detuning and intensity were chosen to give a desired fa 
while maintaining a constant spontaneous emission rate. The range of 0d examined in 
this way was from 0d = 3.3 to 6.6. Taking reflection losses at the cell windows into 
account, over this 0d range the Rabi frequency varied from Q/2n = 34 — 76 MHz with 
corresponding detunings of 645/271 = 315 — 740 MHz, thus 5 » T,Q (where T is the 
linewidth of the atomic transition) was satisfied for all 0d values considered. 

For a chosen 0d the pulse length was held constant (r p = 520 ns), while the pulse 
period was varied from 2.5 /xs to just above the quantum resonance at T « 60 /xs. Thirty 
kicks were delivered to the cloud for each pulse period. The images of the expanded cloud 
were averaged over the dimension perpendicular to the kicking beam to yield momentum 
distributions for the kicked cloud and the mean energy E = (p 2 )/2 was calculated for 
each distribution. High momenta have a large effect on these energy values and as this 
was where the signal dropped for the higher energy kicked clouds, much care was taken 
to reduce the effects of noise. For an experimental run involving a single 0d value, the 
subtracted background was an average from just before and just after the run. Any 
slight fluctuations in background level were accounted for by defining the zero level 
for each momentum distribution via an image taken just before commencement of the 
kicking sequence, omitting the small cloud. The signal-to-noise ratio was on average 
about 100:1 and for each value of T the momentum distribution was measured 5 times. 
To minimise the effects of long term fluctuations in the kicking beam, one fibre end was 
angle cleaved and the beam power was checked and readjusted several times throughout 
a run. These measures reduced the fluctuations in 0d from this source to ~ 1%. 

We now examine the cause of the spread in kicking strengths. Firstly, the finite 
width of the kicking beam results in atoms at different radial positions across the beam 
interacting with laser fields of different intensity. The effect of the spread in kicking 
strengths on the initial Shepelyansky diffusion rates is easily examined and also applies 
to the late-time energies. If p(r) is the 2-dimensional cloud density as a function 
of the radial coordinate r, and 0d( r ) the distribution of kicking strengths, then we 
can approximate the diffusion rate over all of the particles and for a given 0d,max as 
D(fa m T) = D(fa(r),T)p(r)2rrrdr, where D{fa{r),T) is given by Eq. ©. 
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Figure 2. Initial momentum diffusion rates, D', for constant </>d (dash-dotted line) 
and for with a distribution centered on 0d,max = </>d/0.77 with standard deviation 
O.180d (solid line). D 1 is given in experimental units as a function of T for (a) <f)& = 3.3, 
(b) fa — 5.0 and (c) 0d = 6.6. 

Calculating D(0 dmax ,T) for a broad kicking beam with cr c i ou d << 0beam, corresponding 
to constant (pd(r), reproduces the pronounced structure as seen in Fig. [T] and in the 
simulations. But calculating the diffusion rate with the 2:1 beam-to-cloud width ratio 
as used experimentally and with 0d,max = 0d/O.77, gives less accentuated diffusion 
resonance structure which is broader over the pulse period values. Fig. El displays these 
results for a few values of Thus, a spread in d values creates an averaging effect 
which somewhat diminishes the height of the resonance structure and shifts the positions 
of the resonant peaks relative to the case for a single 0d value. 

Secondly, a spread in 0d values is also caused by atoms in different magnetic 
substates of the F = 4 level coupling to the higher energy states with different 
transition strengths, resulting in atoms in different substates experiencing different 
kicking strengths. 

The combination of this effect with that caused by the finite beam width can create 
a spread of 0d values with a standard deviation as large as 20% of the mean 0d value. 
In order to account for this, we performed additional simulations in which the 0d value 
for each trajectory was chosen from a distribution based on the 2:1 beam-to-cloud width 
ratio for our system (assuming that each has a radial Gaussian profile). Also included 
were the relative coupling strengths of atoms in different magnetic substates for the 
experimental detunings 5 = 315,385,485,575,740 MHz for d = 3.3,4.0,5.0,5.9,6.6 
respectively (assuming an equal population of atoms in each substate). 
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Figure 3. Quantum simulations with fixed <j>d (top) and the corresponding 
experimental results (bottom) giving the energy, E', after 30 kicks as a function of 
pulse period with r p = 520 ns, rj = 0.0125 for (a) 0d = 3.3, (b) </y = 4.0, (c) (f>d = 5.0, 
(d) 0d = 5.9 and (e) </>d = 6.6. Additional quantum simulations (middle, solid line) 
take into account the spread in physical kicking strength 0d, as do the analytical 
classical results (middle, dotted line) . The energies arc in experimental units and error 
bars for the experimental and simulation results are shown but are very small. 



5. Experimental Results 

Fig. El shows results from numerical simulations as well as those from experimental 
measurements. The first row shows simulated energies for different, well defined values 
of 0d- All simulations were performed using the Monte Carlo wavefunction method as 
in Ref. ^U]. The simulations are for a system with an initial Gaussian momentum 
distribution of width a p /k = a p /2hki = 4 (this corresponds to a cloud of temperature 
ss 20nK). 

To take into account the unavoidable spread of 0d values in our experiments, 
simulations were also performed in which the kicking strength for each trajectory was 
sampled from the theoretical distribution of 0d experienced by a spherical cloud of atoms 
(see Section^), and the energy for each value of T was taken to be the incoherent average 
of the energies over all such realisations. The results of these simulations are seen in 
the second row of Fig. El Comparison of the first two rows shows that, while the non- 
uniformity of the kicking beam intensity makes the diffusion resonances somewhat less 
distinct, the structure of interest is still clearly visible and thus amenable to experimental 
investigation. 

Comparison of experimental results with those from the simulations (bottom and 
middle rows of Fig. El respectively) shows very good agreement between the two. The 
quantum resonance at T ^ 60 /is is seen to be present for all 0d values. Additionally, 
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while a single broad peak similar to that seen by d'Arcy et al. is found for 0d = 3.3, 
for larger <ftd values more complicated structure is observed. This peak splits into 
two peaks which then diverge from each other, whereupon a third peak rises between 
them. The structure mirrors that in the initial quantum diffusion rate as given by 
Shepelyansky's result in Eq. (jHJ) and shown in Fig. [T] 

There is still some discrepancy in energy values between the experimental results 
and the additional simulations. For = 5.9 and 6.6 the measured energies are, overall, 
lower than expected, while for all 0d they are larger than in the simulations around 
the quantum resonance region of T = 58-68 /xs. The first problem is accounted for by 
realising that in recording the momentum distributions, at some point signals at higher 
momenta fall below the noise level of the CCD. Hence, the measured total energy of 
the cloud is systematically lower than the true energy after kicking. For larger 0d, for 
which energies are higher in general, this problem is particularly pronounced as more 
atoms lie in the wings of the distribution. 

The discrepancy around the quantum resonance currently remains unexplained, 
but could be a systematic effect related to the larger pulse period values in this region. 
Continuing investigations will hopefully resolve this issue. However, we would like to 
emphasise that the energy structure seen at values of T lower than that at which 
quantum resonance occurs is of primary interest to us (experimentally, the quantum 
resonance itself provides a convenient confirmation of the pulse timing calibration). 

Overall, we find very good qualitative agreement between experimental results and 
our simulations. The expected structure is evident and its complexity increases with fa 
as in simulations. In particular, the resonant peaks are seen where they are expected 
for each value of fa. For comparison, the energies of the classical system after 30 kicks 
were also computed by assuming that the diffusion rate for the first two kicks is fa/4:, 
and for subsequent kicks is given by Eq. (JJJ). They were averaged over the same spread 
in 0d values that was used in the additional quantum simulations and are shown in the 
middle row of Fig. El Note that the energies are larger than the quantum simulations, 
as would be expected, and for </>d = 5.9 and 6.6 go off the scale, oscillating around 
E' = 300 and 560 respectively. 

The classical energies clearly exhibit very different behaviour from those measured 
in the quantum system. The difference becomes more marked at higher values of T 
(i.e. when the system behaves more quantum mechanically) where the oscillations in 
the classical energy die down, whilst pronounced diffusion resonances and quantum 
resonances are seen in the quantum system. This can be contrasted with the results of 
Klappauf et al. where, in both the classical and quantum systems, the energy oscillates 
about the quasilinear value and only a relative shift in the peak positions demarcates the 
two regimes. The striking difference between the quantum behaviour and the expected 
classical behaviour in these experiments makes them a thorough testing ground for the 
effects of decoherence on the quantum kicked rotor. Indeed, preliminary investigations 
into the effects of noise on these diffusion resonances [Ej suggest that quantum and 
diffusion resonances are affected differently by noise induced decoherence for a quantum 
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kicked rotor with environmental coupling. 



6. Conclusion 



We have presented experimental and simulation results showing non-trivial behaviour 
in the late-time energy (and thus diffusion rate) of the atom optics kicked rotor 
for values of the pulse period intermediate between those periods associated with 
quantum resonances. The observed resonances are distinct from the quantum resonances 
previously studied in the kicked rotor, being dependent on the kicking strength as well 
as the value of the pulse period. 

The non-trivial features predicted by quantum simulations are matched very closely 
by our experimental measurements. The relationship between the observed structure 
and Shepelyansky's scaling for the quantum correlations is evident. Furthermore, we 
note that the structure observable in the late time energies confirms that in a system 
subject to environmental coupling a dependence of the late-time diffusion rate on the 
initial quantum diffusion rate exists. 

In addition, the observed energy measurements show a very different trend to that 
predicted by classical calculations. With the introduction of noise, these experiments 
should provide a novel testing ground for the efficacy of decoherence in regaining the 
classical limit, and we have begun preliminary studies in this area. 

We thank Maarten Hoogerland for interesting discussions related to our 
experimental investigations. This work was supported by the Royal Society of New 
Zealand Marsden Fund, grant UOA016. 
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